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Abstract 

We show a method to protect quantum states from the disturbance due to the random potential 
by successive rapid manipulations of the quantum states. The quantum states are kept undisturbed 
for a longer time than the case of the simple trapping with a stationary potential. The effective 
potential, which the quantum states feel, becomes uniform when the velocity of the transport is 
sufficiently large. It is also shown that the alternating transport of a Bose-Einstein condensate 
with the driving potential derived by fast-forward scaling theory [Masuda and Nakamura, Proc. 
R. Soc. A 466, 1135 (2010)] can protect it from the disturbance. 



PACS numbers: 37.90.+j, 67.85.d, 81.16.Ta 
I. INTRODUCTION 

Technology to control quantum systems is 
rapidly evolving, and various methods to ma- 
nipulate quantum states have been reported 
in Bose Einstein condensates (BEC)llH4J], in 

n 

quantum computing |5[ and in many other 
fields of applied physics. For many current 
and future technologies the acceleration of 
controls of quantum systems would be impor- 
tant. Methods of the acceleration of quantum 
dynamics and quantum adiabatic dynamics 
or shortcut to adiabaticity have been pro- 
posed, e.g., counterdiabatic protocol js] and 



a, 



frictionless quantum driving 131, invariant- 
based inverse engineering s| and fast-forward 
scaling theory [9l-|lll|. These theories make 



possible to generate target states in short 
time without energy excitations at the final 
time of the manipulations 12] . Recently ap- 
plications of these methods to the controls of 
BEC including the transport have been pro- 
posed theoretically |g[ 13- 181, and been 



demonstrated experimentally 



19 



-|2l|. The 



robustness of these protocols have been in- 



vestigated 



22|. In this paper we show 



a novel application and advantage of the ac- 
celeration of the quantum dynamics. 
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In actual systems there must be noise 

1 



which prevents accurate controls of quantum driving potential for the adiabatic transport 
systems. The influence of the random poten- is reviewed. In Sec JIIII we show the protection 
tial to the static and dynamical properties of of quantum states from the disturbance due 
the Bose-Einstein condensates (BECs) have to the random potential in the large- velocity- 
been studied by using the optical speckle po- limit. In SecJTVl it is numerically exhibited 
tential 



23 



24|. The protection of quantum that the fidelity is kept close to unity by 
states from the influence of the noise and dis- the rapid one-way and alternative transports, 
order would be an important issue as well as The protection of a Bose-Einstein condensate 
rapid controls for accurate manipulations of is also shown numerically, 
quantum states like BECs and for the exten- 
sion of the range of the quantum controls. In jj MODEL 
this paper we use the random potential as the 

noise which deforms the wave function in a We consider a transport of a particle in 

trapping potential neglecting the dissipation one dimension. The driving potential for the 

and decoherence due to the environment un- ideal transport of quantum states without 



ike the quantum decoupling of open systems disturbance was derived, see e.g. [10(. Sup- 
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26|. We show the protection of the quan- pose that ^ (^) is the wave function of an 

turn states from the disturbance due to the energy eigenstate trapped by the stationary 

uncontrollable random potential in the back- potential V (x) in the case without the ran- 

ground by rapid transport of the quantum dom potential. The energy is assumed to be 

states with the use of the fast-forward scaling zero for the simplicity. The potential 

theory. It is shown that the quantum states r} 2 F?(t) 

V(x, t) = V (x- R{t)) - ^mx. (1) 
are kept undisturbed for a longer time than at 2 

the case of the simple trapping with a station- can translates the quantum state without en- 

ary potential because the effective potential, ergy excitation at the final time of the manip . 

which the quantum states feel, becomes uni- uktion Q The ch&nge ^ R(t) ig the dig _ 

form when the velocity of the transport is suf- placement of the trapping potential and the 

ficiently large. It is numerically exhibited in wave function . The first term in Eq & cor . 

one dimension that the alternating transport responds to the translation of the trapping 

of a Bose-Einstein condensate can protect it poten tial. The second term is the additional 

from the disturbance. potential which is spatially linear. The wave 

In SecJTT]we represent the model. And the function of the transported state ^ff is rep- 



x exp 



resented as the random potential in the case of the con- 

stant velocity, R — v, in the large-velocity- 

V FF (x,t) = V (x-R(t)) 

limit. In the analysis we use the moving 

:(*)— x (2) . . . . . . . 

h J frame which accompanies with the trapping 

where R denotes the time-derivative of R. potential. In the moving frame the third term 

R(t) is the velocity of the translation, in Eq.flU) vanishes and the Hamiltonian is 

%> FF (x,t) is a solution of the Schrodinger represented by 
equation: 2 

„»„ ,„ „ „ ^ = £ + V. W+ V,(x + rt ). (5) 

The trapping potential is stationary in the 

The additional phase in the wave function in . 

moving frame while the random potential is 

Eq.((2D vanishes everywhere when the quan- . . . TTr 

moving with the constant velocity. We ex- 

tum state is stopped and i? = 0. , , T , , 

pand the state V by the energy eigenstates 

Now we consider the translation under the . . „ , TT .„ TT 9 

\j > of the Hamiltonian: H = p /2m + 

random potential VJx). We assume that the Tr , , . , 

Vo(x). Ihe state is represented as 

random potential is time-independent. In 

general the random potential can cause the >= ^^ a j(t)\j >> (6) 

disturbance of the wave function. The to- 3 

tal Hamiltonian with the driving potential is where H \j >= Ej\j > and < j\k >= 5j k . 

represented as Ej is the energy of the jth state. \j > is 

v 2 , . time-independent. The Schrodinger equation 

H = f- + V (x-R{t)) 

^ m leads to the equations of the coefficients: 

^ LmX + V r( X )- ( 4 ) 

d taj (t) = ~E jaj (t) - % T^V jk {t)a k {t), (7) 
In the following sections it is shown that the k 

disturbance due to the random potential is , T/ ,.\ ■ , n 11 

^ where Vj k {t) is denned by 

restrained by the rapid transport. 

V jk (t) = < j\V r (x + vt)\k > 

/oo 
(p*(x)(p k (x)V r (x + vt)dx,(S) 
-00 

ITY LIMIT 

with 4>j(x) =< x\j >. We suppose that the 
We show that the fast transport of the initial state is the nth energy eigenstate \n >, 
quantum states can reduce the influence of that is, %(0) — Sj n . dj(t) satisfies the integral 



equation: where u mn = (E m — E n )/H. For example the 



aj (t) = e-&i% 



4E f ds e-^ {t - s) V 3k (s) ai 
n k J o 

(9) 

For m 7^ n we obtain the 1st order ap- 



gth integral with respect to s of the last line 
in Eq. (TTUI) is represented as 

pqAt 

e ^ mn qAt / dsV mn { S ) 
J(q-l)At 

pqAt poo 

= e ^ mnqAt / ds dx 

J(q-l)At J -oo 

proximation of a m (t) by substituting j = X(j)* m {x)(j) n {x)V r {x + vs) (13) 



m and the Oth order solution dfc(s) 
5kn exp[— iE k s/h] in Eq.((9]) as 



a m (t) ~ j/ s e~^- s) pAt 



The integration with respect to s is rewritten 
with r = x + v s as 



-i e -**»* = At iA7 / K(r)rfr - (14) 



X 



/ ds e^ Em ~ En ^ s V mn (s). (10) We define the effective potential K(^,<7) by 



We assume that a finite number of the energy K(x, q) = / V r (r)dT, (15) 

vAt J x+v ( Q -i)At 

eigenstates can have the dominant contribu- 



tion to the transition from the nth eigenstate 
in the 1st order approximation and the oth- 
ers are negligible because of the sufficiently 
small V mn . Hereafter we focus on the en- 
ergy eigenstates which may have the domi- 
nant contribution. We divide the time inte- v ^ > ^ 

we have 



x+vqAt 
x+v(q-l)At 

which is the average of the random potential 
in the interval: Ax = vAt. Suppose that / 
is the length of vAt with which V r (x) can be 
regarded as uniform in the region that wave 
function is localized. Thus if 



gral in Eq. fjlQ]) into the intervals as 

t rAt p2At ft 

ds = ds + / ds H — • + / ds 



iAt 

(11) ~ mw - ~T 



Jt-At a m (t) ~ c ~^ Emt c iuJm "i At 

q 



/oo 
dx (j)* m (x)(j) n (x)V r (x,q) ~ 0,(17) 
-oo 

E n )s/h] can be regarded as constant in the 

due to the orthogonality of <j) n and <j) m . 

interval for any m, that is, 

Therefore in the large v limit, we see no tran- 
u mn At <C 1, (12) sition among the energy eigenstates. From 



the conditions in Eqs.f ll2p and ([16]) We ob- due to the random potential, while they coin- 
tain a criterion of v: cide with each other at the initial time. The 

random potential has the rectangular form 

v > u mn l (18) 

with the width rj as shown in FigHJ The ran- 
with which the level transitions due to the dom potential takes the value from -W/2 
random potential do not occur. to w / 2 . We drive the wave function in x- 

I 

We numerically exhibit the protection of 
quantum states by the one-way translation 
with the constant velocity and the alternat- 
ing translation in one dimension. 

0.025 




A. Translation with constant velocity 

We numerically simulate the transport 
with the constant velocity under the random 
potential. The fidelity of the quantum state 
is calculated during the time-evolution. We 
chose the trapping potential Vq in Eq.(j4]) as 
the harmonic potential: 



V (x) 



mu 



-x 



(19) 



The initial state is taken as the ground state 
in the harmonic potential multiplied by a 
phase factor as 



tf(z,0) 



mu 



exp 



mu 



2h 



-x 



■ m 
iR—x 
n J 



(20) 



where R is the velocity of the transport and 
is constant. ^>(x,t) is deformed from the ex- 
actly transported state ^FF{x,t) in Eq.(J2]) 



FIG. 1: Snapshots of the time-evolution of l^j 2 
under random potential for R = 0. Red and 
blue lines correspond to j^j 2 and the random po- 
tential, respectively. The parameters are taken 
as m = 1, h = 1, u = 256, 7] = 0.065 and 
W = 1000. 

direction by moving the harmonic trap with 
the constant velocity (see Eq.fllJ)). In the nu- 
merical simulation we use the moving frame 
which accompanies with the trapping poten- 
tial. The corresponding Hamiltonian is given 
by Eq.©. 

The snapshots of the time-evolution of 
|\1/| 2 under random potential are shown in 
FigUJ for R = which corresponds to the 
simple trapping with the stationary poten- 
tial. The harmonic potential is not shown 



in the figure. The parameters are taken as 
m — 1, H — 1, oj — 256, 7] = 0.065 and 
W = 1000. The wave function is deformed 
due to the random potential. The time de- 
pendence of the fidelity is shown in FigfS] for 
various values of R. The fidelity is defined by 
| < ^ FF (t)\^(t) > | where \^ FF > in Eq.© 
is the exactly transported state without the 
disturbance of the random potential. The fi- 
delity is averaged over the dynamics with 100 
different random potentials. We see the de- 
crease of the fidelity with time due to the dis- 
turbance by the random potential for R = 0. 
The fidelity for R = 10 and 30 is lower than 
that of R = 0. This is because that the rel- 
ative time- dependence of the random poten- 
tial enhances the energy transition. However, 
for the sufficiently large velocity, R > 100, 
the decrease of the fidelity is apparently re- 
strained compared to the case of R = 0. 
The continuous transport reduces the influ- 
ence of the random potential as theoretically 
predicted in the previous section. The vari- 
ance of the fidelity is about 0.01 for R = 100. 
The variance tends to be smaller for larger R. 

V r in Eq. ffTS]) is regarded as the effective 
potential that the quantum state feels. To 
show the property of the effective potential 
we calculate 

V eff (x,5l) = - V r (r)dr, (21) 




FIG. 2: Time dependence of fidelity in one-way 
driving. The fidelity is averaged over the dy- 
namics with 100 different random potentials. R 
is the velocity of the centre of trapping potential. 
Other parameters are the same as FigfJJ 

which corresponds to V r for q = 0. The 
distance 51 corresponds to v At. The effec- 
tive potential V e ff(x,Sl) for the various val- 
ues of 51 are shown in FigfSl For large 51 
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FIG. 3: Effective potential V e ff(x,5l) for 51 = 
0,0.1,1,10. V e ff for 51 > are shifted upward 
for the comparison. Other parameters are the 
same as Fig(TJ 



(51 = 1 and 10) the effective potential be- accompanies with the original trapping po- 

come smooth and uniform compared to the tential Vq. 

original random potential. This property of The time dependence of the fidelity is 

the effective potential explains the decrease shown in FigHfor various values of u R . The 

of the influence of the random potential. The corresponding value of the time-average of 

distance 51 = 1 corresponds to R = 256 if \R\ denoted by < \R\ > are shown in the fig- 

we relate the distance 51, the time-interval U re. We put L = 2. Other parameters are the 

At and the frequency u by 51 = RAt and same as FigJU The fidelity is averaged over 

At = 1/uj. the dynamics with 100 different random po- 
tentials. The variance of the fidelity is about 

B. Alternating translation °- 015 for < 1^1 >= 4000 - The variance tends 

to be smaller for larger < \R\ >. It is seen 

In actual systems it is impossible to keep that the alternating translation with the suf- 

translating the quantum state in one direc- ficiently large frequency can reduce the influ- 

tion. Here we show that the alternating ence Q f random potential. The rapid falls of 

translation also protects quantum states from the fidelity occur at the time when the trap- 

the disturbance due to the random potential ping potential is turned and starts to move in 

in one dimension. The trapping potential the opposite direction. Since the velocity be- 

V (x) is the harmonic potential as Eq.JTSD, come s small at such time, the quantum state 

and the initial state is the ground state given is affected by the random potential and the 

by Eq.(f20D with R = 0. We choose the time fidelity decreases, 
dependence of R(t) as 

R(t) = —{1 — cos(ojRt)}. (22) C. Bose-Einstein condensates 

We continuously translate the quantum state Here we apply the method to reduce the 

back and forth by translating the trapping influence of the random potential to Bose- 

potential and simultaneously tuning the spa- Einstein condensates (BECs). We assume 

tially linear potential as Eq.flTJ). The center that the system is governed by the Gross- 

of the wave packet is oscillated between x = Pftaevskii (GP) equation: 
and x = L periodically. The random poten- 

& d 2 ^ w lTl2T 

tial is the same as Sec flVAl In the numerical ^— = + V(z, + c|W| , 

simulation we use the moving frame which (23) 
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FIG. 4: Time dependence of fidelity in the alter- 
nating translation for various values of < \R\ >. 
The fidelity is averaged over the dynamics with 
100 different random potentials. Other parame- 
ters are the same as FigHJ 

where ^>(x,t) is the macroscopic wave func- 
tion, and c is the coupling parameter. It has 
been shown that the BECs are also trans- 
ported without energy excitation by the same 
driving potential in Eq. (JT]) [lOj, |l8|, |27J . Let us 
suppose that a BEC wave packet is trapped 
by the harmonic potential subjected to the 
random potential which is the same as the 
previous sections. The initial state is the 
ground state obtained numerically in the sta- 
tionary harmonic potential without the ran- 
dom potential. We alternately translate the 
wave packet by using the driving potential. 
The time dependence of R(t) and the pa- 
rameters are same as Sec JIV Bl In Figj5] the 
time dependence of the fidelity is shown for 
c = 10. The fidelity shows the similar time- 



dependence to the case with c = 0. The re- 
sults show that the disturbance of BECs can 
be suppressed by the repeated driving of the 
wave packet, 
l 
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FIG. 5: Time dependence of fidelity in alternat- 
ing translation for c = 10. Other parameters are 
the same as Figf4j 

V. SUMMARY AND DISCUSSION 

We have presented the protection of quan- 
tum states from the disturbance due to the 
random potential by the continuous trans- 
port with the use of the driving potential 
derived by the theory of the acceleration of 
adiabatic dynamics. We analytically showed 
that the fast translation of the quantum 
states can suppress the influence of the ran- 
dom potential because the effective potential 
which the quantum states feel becomes uni- 
form when the velocity of the transport is suf- 
ficiently large. We emphasize that the veloc- 
ity of the transport does not have to be fast 



S 



enough so that the effective potential van- consider the dissipation and decoherence ef- 

ishes. fects. Applications of the present method in 

We have numerically exhibited the protec- the system with the dissipation and decoher- 

tion of the quantum states by the continuous ence would be studied in the future also with 

translation with the constant velocity and the effect of the time- dependence of the ran- 

the alternating translation of the the ground dom potential. We assumed that the driving 

state in the harmonic potential. The decrease potential is controlled without error, while for 

of the fidelity is clearly restrained by the fast- rapid manipulation the control of the poten- 

driving compared to the simple trapping with tial itself can cause additional noises in actual 

a stationary potential, while the protection systems. The investigation of the robustness 

effect is not monotonously increased with the of the present technique is a future issue, 
velocity in the small-velocity-region. It has 
been shown that the same technique is effec- 
tive also for the protection of BECs from the Acknowledgments 
influence of the random potential. 
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